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AUXILIARY LOOP FOR REMOVAL OF DECAY HEAT FROM 
A SPACE POWER REACTOR 
by George E. Turney and Edward J. Petrik 
Lewis Research Center 

SUMMARY 

An auxiliary loop using sodium as a coolant was investigated as a means for remov- 
ing decay heat from a proposed 2. 17-megawatt space power reactor. This loop consists 
of an auxiliary radiator, an auxiliary heat exchanger, a liquid sodium pump, and the 
necessary flow control valves and piping. 

Based on a selected set of operating conditions and assumptions, minimum areas 
were determined for the auxiliary loop heat exchanger and radiator. The minimum 

auxiliary heat exchanger surface area required to transfer the decay heat from the re- 

2 2 

actor loop to the auxiliary loop was 0. 64 ft (0. 0594 m ). The minimum auxiliary radi- 

2 2 

ator surface area required to reject the decay heat was 37. 7 ft (3. 50 m ). Based on 
these areas, the maximum temperature of the auxiliary radiator was 1760° R (978 K). 

The configurations and minimum areas for the auxiliary loop heat exchanger and 
radiator were used to determine sodium flow rates and temperatures throughout the 
auxiliary loop as a function of time after reactor shutdown. The sodium flow rates and 
temperatures were determined for preshutdown reactor operating times of 1 year and 
1 day. 

It was concluded from this preliminary study that the size (and/or mass) of the 
auxiliary loop is relatively small. 


INTRODUCTION 

In the coming years, the electric power requirements for our nation's projected 
space programs will continue to increase. To meet these projected requirements, the 
Lewis Research Center has been participating in a technology program aimed at a rela- 
tively high-powered nuclear Brayton space powerplant. 

The heat source being considered for the Brayton cycle space powerplant is a com- 



pact, fast- spectrum nuclear reactor. The design thermal output of the reactor is 2. 17 
megawatts, and the design operating lifetime for the reactor and powerplant has been set 
at 50 000 hours. 

An area of concern in the overa.ll design of this powerplant is the removal of decay 
heat from the nuclear reactor following a system shutdown Although shutdowns are un- 
desirable, particularly in a space power system, they must be considered in the system 
design. And provisions must be made for removing the reactor decay heat following a 
shutdown. 

In this report, an auxiliary radiator loop with sodium coolant is investigated as a 
means for removing reactor decay heat. The objectives of this preliminary investigation 
are (1) to determine the size of an auxiliary radiator which will reject the decay heat and 
operate at an acceptable temperature level during the shutdown period, (2) to determine 
the size of an auxiliary heat exchanger which will transfer the decay heat from the re- 
actor primary loop to the auxiliary radiator loop, and (3) to evaluate the sodium flow 
rates and temperatures in the auxiliary radiator loop as a function of time after shutdown. 

An important consideration related to shutdown is the ability to restart the system 
and return it to the design operating level in a minimum amount of time. Based on this 
consideration, a mode of operation was chosen in which the reactor coolant temperatures 
were maintained at their respective design values during the period of shutdown. 

DESCRIPTION OF POWERPLANT 

A simplified diagram of a conceptual nuclear Brayton space powerplant is shown in 
figure 1. This system, in its present design configuration, has a single liquid-metal 
(lithium) primary loop, one or more complete inert-gas power conversion loops, and a 
main radiator loop for waste heat rejection. 

The path of coolant flow in each of the three main loops is indicated in figure 1. 
During normal operation, the heat gained by the liquid lithium as it flows through the 
reactor is transferred to the gas power conversion loop by a heat- source heat exchanger. 
The heated inert gas is then expanded through a turbine. The work produced by the tur- 
bine drives a compressor and an electrical alternator. The exhaust from the turbine is 
cooled as it flows, first through the recuperator, and then through the waste heat ex- 
changer . 

Figure 2 shows a cutaway view of the reactor core design for this space power sys- 
tem. A detailed description of this reactor design is given in reference 1. The system 
design point operating conditions which are relevant to this study are listed in table I. 
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Gas power conversion loop 



Figure 1. - Schematic diagram of a three-loop, recuperated, nuclear Brayton system. 
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Figure 2. - Nuclear Brayton space power reactor. 


TABLE I. - DESIGN POINT OPERATING CONDITIONS 


Reactor thermal power, mW 
Lithium temperature at reactor inlet, ° 
Lithium temperature at reactor outlet. 
Lithium flow rate, Ibm/sec (kg/sec) 


2. 17 

R (K) 2100 (1167) 

°R (K) 2200 (1222) 

20.7 (9.39) 




AUXILIARY LOOP 


General Considerations 

Following a shutdown of the nuclear reactor, decay heat continues to be generated in 
the reactor core at a rate which depends primarily on (1) the reactor's initial operating 
power level, and (2) the period of time during which the reactor was in operation before 
the shutdown. In a space power reactor, the decay heat must ultimately be rejected to 
space by some type of waste-heat radiator. 

It may be both impractical and undesirable to use the main radiator loop of the sys- 
tem (fig. 1) as a primary scheme for rejection of decay heat. It is conceivable that a 
shutdown may be caused by a failure in the main radiator loop. In such case, it may be 
impossible to use the main radiator . Furthermore, it may be desirable in the system 
design to have a method for rejection of decay heat which is independent of the main sys- 
tem operation . 


Auxiliary Loop Description 

One scheme for rejecting the decay heat is to use a liquid- metal auxiliary radiator 
loop. Figure 3 shows a sodium- cooled auxiliary loop which is connected to the reactor 
primary loop by an auxiliary heat exchanger . The auxiliary loop consists of a high- 
temperature auxiliary radiator , an auxiliary heat exchanger, a liquid- sodium pump, and 
the necessary flow control valves and piping. Sodium was selected as the coolant in the 
auxiliary loop because it has a reasonably low melting point and also a relatively low 
vapor pressure for the temperatures which we will be considering. 

During normal system operation, nearly all the heat generated in the reactor is 
transferred to the gas power conversion loop. But when a shutdown occurs, the decay 
heat generated in the reactor must be transferred to the auxiliary loop. To accomplish 
this, the flow of lithium in the primary loop must be bypassed around the heat- source 
heat exchanger; and the sodium flow rate in the auxiliary loop must be increased to 
transfer the decay heat. Hence, the auxiliary loop undergoes a transient startup imme- 
diately following the reactor shutdown. 


Operating Assumptions 

To investigate the operation of this loop, some assumptions were made regarding 
the initial conditions in the loop and the mode of operation of this loop during the reactor 
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Auxiliary loop 
(sodium coolant) 



Figure 3. - Primary loop with auxiliary radiator loop for decay heat rejection. 
(See appendix for number designations.) 


shutdown transient. For the preliminary investigation presented herein, we assume that 
the shutdown is not the result of a serious failure and that, as a consequence, the sys- 
tem is capable of being restarted and returned to its design operating condition. 

As a point of interest, the restarting of this fast-spectrum reactor is, in some re- 
spects, less complicated than that of a thermal- spectrum reactor. This is because the 
xenon- 135 and samarium- 149 concentrations, which build up after shutdowns, do not act 
as strong absorbers (i. e. , their absorption cross sections are small) for high-energy 
neutrons. And the problem of overriding the poisoning effects of these isotopes is far 
less serious in this reactor. 

The major assumptions used in this preliminary analysis and related to the operat- 
ing conditions are listed below: 
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( 1 ) The lithium temperature at the reactor inlet is maintained constant at 2100° R 
(1167 K) during the shutdown. (This is the design point operating temperature (table I). ) 

(2) The lithium temperature at the reactor outlet is held constant at 2200° R (1222 K) 
during shutdown. (This is also the design point operating temperature (table I). ) 

(3) The initial temperature of sodium in the auxiliary loop at shutdown is 100° R 
(55. 6 K) above the melting temperature of sodium, or 770° R (428 K). 

(4) In order to keep the auxiliary loop temperatures reasonably uniform after shut- 
down, we assume the mixed-mean temperature difference across the auxiliary heat ex- 
changer to be constant at 100° R (55. 6 K); that is, (T^^^ g - Tj^^ 5) = 100° R (55. 6 K). 
(All symbols are defined in appendix A. ) 

(5) A maximum allowable sodium temperature of 1760° R (978 K) at the auxiliary 
radiator inlet is assumed. 

Assumptions (1) and (2) are consistent with a requirement that the system be capable 
of being restarted in a minimum time. Maintaining the primary loop temperatures at 
their respective design values will minimize the time required to return the system to 
its design operating level. 

The selection of an auxiliary loop temperature of 770° R (428 K) (assumption (3)) was 
based on several considerations. First, this initial (or standby) temperature level of 
770° R (428 K) gives a reasonable margin of safety against freezing of sodium in the 
auxiliary loop. At the same time, this temperature is low enough during standby oper- 
ation that the heat rejected by the auxiliary radiator (which is a parasitic thermal load) 
is small relative to the reactor design thermal power. 

Assumptions (4) and (5) limit the operating temperature and the differential tempera- 
ture in the auxiliary loop. These values should not be considered as absolute limits. 
They were selected somewhat arbitrarily for this preliminary investigation in order to 
place some constraint on the operating conditions in the auxiliary loop during the decay 
heat removal period. 


Auxiliary Loop Flow Models 

In order to make a preliminary investigation of this loop, simplified flow models 
were developed to represent the auxiliary radiator and auxiliary heat exchanger. The 
auxiliary loop material (including the auxiliary heat exchanger, auxiliary radiator, and 
piping) is assumed to be tantalum. All tantalum physical properties used were taken 
from reference 2. A description of the flow models is given below. 

Auxiliary heat exchanger . - The auxiliary heat exchanger is a counterflow type, 
consisting of two concentric cylindrical ducts. Figure 4 is a sketch of the auxiliary heat 
exchanger. Lithium is assumed to flow through the inner passage, and sodium flows 
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counter currently in the annular passage. Based on preliminary estimates of the area 
requirements of this heat exchanger, a cross-sectional configuration was selected. The 
cross-sectional configuration is described in table II. 

Auxiliary radiator . - The auxiliary radiator is made up of a number of coplanar 
parallel tubes with central connecting fins. For this study, the effectiveness of the fins 
was assumed to be 100 percent. A sketch of the auxiliary radiator configuration is 
shown in figure 5. 

Preliminary estimates were made to determine the area requirements and a reason- 
able planform configuration for the auxiliary radiator. From these estimates, a cross- 
sectional configuration was selected. It is based on the assumption that heat is rejected 
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Number of tubes 

Outside diameter of tubes, in. (cm) 

Inside diameter of tubes, in. (cm) 

Spacing between tubes (or central fin width), in. (cm) 
Central fin thickness, in. (cm) 


24 

1 (2.54) 
0.875 (2.22) 
1.055 (2.68) 
1/16 (0. 159) 


to space from both sides of the auxiliary radiator. A description of the auxiliary radi- 
ator cross-sectional configuration is given in table III. 

The central fin width (1. 055 in. or 2.68 cm) was chosen for convenience. With this 
fin width, the surface area of all tubes is equivalent to the product of total radiator sur- 
face area (which includes tubes and connecting fins) and the configuration factor from the 
auxiliary radiator to space. 


ANALYSIS 

Heat -Transfer Correlations 

The correlations used to compute the rate of heat flow in the auxiliary loop are de- 
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scribed below. The heat- transfer coefficient on the lithium side of the auxiliary heat ex- 
changer was computed from the correlation (ref. 3) 

N“Li = »26(Re 8 (1) 

The lithium properties used in this study were obtained from reference 4. 

The sodium heat-transfer coefficient of the annular flow passage was determined 
from the equation (ref. 5). 


Nu 


Na 



5.H-0:028(pe^J'’-8 


( 2 ) 


The physical properties of sodium were taken from reference 6. 

The overall heat-transfer coefficient for the auxiliary heat exchanger, based on the 
surface area of the sodium side, is given by 


U, 


/D. 


D, 


D. In - 


b 


b 




°a VhLi^ 


2k. 


m 


m 


(3) 


In the auxiliary radiator, an infinite heat-transfer coefficient was assumed on the 
liquid- sodium side of the flow passages. This assumption is discussed later in this 
section. 


Analytical Equations 

Primary loop . - The lithium flow rate required to maintain a constant temperature 
difference of 100° R (55. 6 K) across the core at any time t after shutdown is given by 

w, . = § (4) 

S,Li<TLi,l-TLl,2) 

where T^^ ^ = 2200° R (1222 K) and Tj^. ^ = 2100° R (1167 K). In equation (4), Q is 
the decay heat generation rate, which varies with time after shutdown. 

Auxiliary loop. - The total flow rate of sodium needed to maintain a constant mixed- 
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mean temperature difference of 100° R (55. 6 K) in the auxiliary loop at any time t after 
shutdown is given by 


^Na, 5 “ ^Na, 8 ^ 




p,Na^%a, 8 " %a, 5^ 


( 5 ) 


From operating assumption number (4) in the section AUXILIARY LOOP, 

(Tnu 8 ” '^Na 5^ constant and equal to 100° R (55. 6 K) throughout the shutdown period. 

The decay heat generation rate and the overall lithium - to- s odium temperature differ- 
ence across the auxiliary heat exchanger change with time following a shutdown. Conse- 
quently, the flow rate of sodium in the auxiliary heat exchanger Wj^^^ rj must also change 
with time. The sodium flow rate needed to transfer the decay heat from primary loop to 
auxiliary loop at any time t after shutdown was determined from the equation 


w- 


Q 


EC .T (Tt • 1 - T.t c) 
p, Na'' Li, 1 Na, 5^^ 


( 6 ) 


The sodium flow rate which bypasses the auxiliary heat exchanger is the difference be- 
tween the total sodium flow rate and the sodium flow rate which goes through the auxil- 
iary heat exchanger; that is, 


^Na, 6 - ■^Na, 5 " ^Na, 7 


(V 


The effectiveness of the counterflow auxiliary heat exchanger (E in eq. (6)) is also 
time variant, and is defined by the following expression (ref. 7): 


1 - exp 


E =• 



L ^pNa^Na, 7\ 

_S,Na^Na, 7 

' '^pLi*LI ) 


^pNa^Na, 7 ^ 

CpLi^Li > 


exp 


Ub\ 


^pNa^Na, 7 
CpLi^Li /J 


^pNa'^Na, 7 


( 8 ) 


The sodium temperature at the outlet of the auxiliary heat exchanger was determined 
from the equation 


"^Na, 7 ~ "^Na, 6 


Q 


^Na, 7^pNa 


(9) 


For calculational purposes, the flow passages in the auxiliary radiator were con- 
sidered to be made up of a number of equal length segments arranged in series. (The 
actual cross-sectional configuration of this radiator is described in table III. ) A sketch 
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Figure 6. - Analytical model used to represent auxiliary radiator. 


of the analytical model for the auxiliary radiator is shown in figure 6. As indicated in 
the figure, each tube in the auxiliary radiator was assumed to be made up of "n" seg- 
ments. A segment length of 1. 8 inches (4. 57 cm) was used in the analysis. The reason 
for using a segment length of 1. 8 inches or 4. 57 cm is given later in this section. 

The following assumptions were used in the auxiliary radiator calculations: 

(1) The temperature drop through the auxiliary radiator tube wall is negligible. 

(2) Heat conduction along the length of the auxiliary radiator is negligible. 

(3) Only the thermal capacitance of the auxiliary radiator (which includes tubes and 
fins) is considered. 

(4) Sink temperature is constant and equal to 400° R (222 K). 

(5) Radiator surface emissivity e is equal to 0. 80. 

With steady flow in the auxiliary radiator, an energy balance on a typical segment 
(e.g. , segment j) is 


''Rate of change of\ 
sodium internal \ 
energy across j 
segment j / 


'Rate of heat trans- 
fer from sodium 
to radiator wall 
in segment j 



Rate of change in'' 
sensible heat of 
radiator seg- 
ment j 


( Rate of heat re- 
jection by radi- 
ator segment j 
to sink 
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Or, in equation form, 


■w 


C HT 
Na, 4 pNa Na 


i,j “^R 


dT 

+ ae AA„fT 




dr 


o\ Na, i 


( 11 ) 


In equation (11) the subscripts i and o on AA refer to the inside and outside surfaces 
of the radiator tubes (see fig. 6). 

With liquid- sodium flow, the wall-to-bulk temperature difference in the auxiliary 
radiator tubes is relatively small. (For a radiator temperature of 1700° R (944 K), the 
wall-to-bulk AT with turbulent flow is less than 10° R (5. 56 K). ) Therefore, we can 
neglect the liquid-metal film resistance and assume that Tj^ ^ = Tj^^^ j without introduc- 
ing a significant error. And equation (11), in finite -difference form, becomes 


^Na, 4SNa^%a, j,in 


'Na 


,) = AM.C 
,],ouR ] 


AT 


M 


Na, j 


+ ere 


aA^ 


At 


( rp4 

rNa,j 


( 12 ) 


where 



T + T 

Na, j,in Na, j , out 


2 


(13) 


The sodium flow rate in the bypass around the auxiliary radiator is given by 


^Na, 3 ~ "^Na, 5 


w 


Na, 4 


(14) 


The mass of each radiator segment (i. e. , the value of AM in eq. (11)) was deter- 
mined from the cross-sectional configuration listed in table III along with the assumed 
length of the segments. As stated earlier in this section, only the mass of the auxiliary 
radiator (which includes tubes and fins) was considered in the study; and the thermal 
capacitance of the auxiliary heat exchanger, inlet plenums, valves, and piping were 
ignored. 

The time interval At used in equation (12) is equal to the auxiliary loop transport 
delay (or delay time). The transport delay is defined as the time required for the sodium 
coolant to traverse the auxiliary loop. Expressed as an equation, 

_ Total volume of auxilia.ry loop 
Sodium volumetric flow rate 
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or 


^ (Volume of auxiliary radiator) + (Volume of auxiliary loop flow lines) 

^Na, 5 

^Na 

q 

A volume of 0. 260 cubic feet (0. 02835 m ) was assumed for the auxiliary loop flow 
lines. This volume is equivalent to that of a 1. 5-inch- (3. 81-cm-) inside-diameter pipe 
with a length of about 21 feet (5.1m). 

The volume of the auxiliary radiator is equal to the product of the auxiliary radiator 
flow area and length. The flow area is obtained from table III, and the procedure used to 
determine the length of the auxiliary radiator is described later in this section. 


Analytical Procedure 

A digital computer code was written to analyze the decay heat removal from the pri- 
mary loop of the nuclear space powerplant. For purposes of analysis, the removal of 
decay heat from the primary loop can be considered to take place in two sequential steps. 
First, for some period of time immediately following a shutdown, most of the decay 
power goes to increasing the sensible heat of the sodium- cooled auxiliary loop. During 
this period, the sodium flow rates through and around the auxiliary heat exchanger 
(fig. 3) are regulated such that the mixed-mean sodium temperature rise across the 
auxiliary heat exchanger is 100° R (55. 6 K). (This is consistent with operating assump- 
tion (4) in the section AUXILIARY LOOP. ) Throughout this auxiliary loop heat-up period, 
the total flow of sodium passes through the auxiliary radiator. And the auxiliary loop 
temperature is raised from its initial value of 770° R (428 K) to some selected final oper- 
ating temperature. (In the section AUXILIARY LOOP, a maximum allowable sodium tem- 
perature of 1760° R (978 K) at the auxiliary radiator inlet was assumed. ) 

The second part of the decay heat removal period begins when the sodium tempera- 
ture at the auxiliary radiator inlet reaches its final operating temperature. Then, the 
mixed-mean sodium temperatures in the auxiliary loop (Tj^^ g and Tj^^ g in fig. 3) 
are held at their respective levels by regulating the flow rates of sodium through and 
around the auxiliary heat exchanger and auxiliary radiator. A sketch of the typical his- 
tory of the sodium temperature at the auxiliary heat exchanger inlet following a shutdown 
is shown in figure 7. 

In the following paragraphs, we describe the procedure used in the decay heat re- 
moval analysis. Decay heat generation rates for the space power reactor were calculated 
for several full- power operating time periods . The method of calculating the decay heat 
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Figure 7. - Typical history of sodium temperature at auxiliary heat exchanger 
inlet as function of time after shutdown. 


rates is described in reference 8. From these calculated values, an upper and lower 
limit for decay heat generation rate was determined. The upper limit corresponds to a 
reactor full-power operating period of 1 year; the lower limit corresponds to a full- 
power operating period of 1 day. 

An iterative technique was used to determine the minimum areas required by the 
auxiliary heat exchanger and auxiliary radiator. The upper limit of decay heat genera- 
tion was used as the basis for these calculations. The calculations were made by picking 
trial values of area (or length) for the auxiliary loop heat exchanger and radiator. 

With these trial values, flow rates and temperatures were calculated in the auxiliary 
loop as a function of time after shutdown. 

The temperatures and flow rates obtained were then used as guidelines to update the 
areas (or lengths) of the auxiliary heat exchanger and auxiliary radiator. The calcula- 
tions were then repeated with the updated areas. This calculational procedure (in which 
successive computations were made with updated values of area or length) was continued 
until the final areas (or lengths) obtained were minimum values consistent with the estab- 
lished temperature limitations for the auxiliary loop. 

The flow rate in the auxiliary heat exchanger Wj^g^ rj was calculated as a function of 
time after shutdown from equations (6) and (8). An iterative technique was used to solve 
these equations for each trial value of auxiliary heat exchanger area. 

The sodium temperature at the auxiliary radiator outlet Tj^g ^ was determined by 
numerically integrating equation (11) from auxiliary radiator inlet (segment 1) to outlet 
(segment n). As stated earlier in this section, the total sodium flow rate passes through 
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the auxiliary radiator during the auxilia,ry loop heatup period. But when the sodium tem- 
perature at the auxiliary radiator inlet reaches 1760° R (978 K), the decay heat genera- 
tion rate must be balanced by the net rate of heat lost by the sodium in the auxiliary radi- 
ator. This criterion was used to determine the minimum area (or length) of the auxiliary 

radiator . 

After the sodium temperature at the auxiliary radiator inlet reaches 1760° R (978 K), 
a part of the total sodium flow must be bypassed around the auxiliary radiator in order to 
maintain the proper heat rejection rates . The flow rates through and around the auxil- 
iary radiator were calculated from equations (12) and (14). An iterative procedure was 
used to determine these flow rates as a function of time. 

The numerical open- loop integration of equation (12) is subject to a cumulative error. 
The significance of the error depends primarily on (1) the size and/or number of seg- 
ments used to represent the auxiliary radiator, and (2) the number of numerical integra- 
tions carried out. 

The magnitude of this cumulative error was investigated by using radiator segments 
with lengths ranging from 0. 72 to 7.2 inches (1. 83 to 18. 3 cm). From this investigation, 
it was determined that a segment length of 1. 8 inches (4. 57 cm) gives a reasonably accu- 
rate representation of the auxiliary loop variables, provided the shutdown times used in 
the analysis are relatively short (i.e. , on the order of 3 hr (10 800 sec) or less). 

All calculations required in this quasi- steady- state analysis were performed on a 
high-speed digital computer . In table IV, we have summarized the values used for the 
input variables in this analysis. These input values have been stated and discussed in 
previous sections of this report. 


TABLE IV. - SUMMARY OF VALUES USED IN ANALYSIS 


Radiator surface emissivity, e 

0. 80 

Sink temperature, T^, °R (K) 

400 (222) 

Lithium temperature at reactor inlet, °R (K) 

2100 (1167) 

Lithium temperature at reactor outlet, Tj^. g, °R (K) 
Initial temperature of auxiliary loop prior to 

2200 (1222) 

shutdown, °R (K) 

Sodium mixed-mean temperature difference, 

770 (428) 

%a,8-%a,5’ ® 

Maximum allowable sodium temperature at 

100 (55.6) 

auxiliary radiator inlet, g, °R (K) 

1760 (978) 
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RESULTS AND DISCUSSION 


Decay heat generation rates in the nuclear space power reactor were calculated as a 

c; X ci.Li.x_; x_;i xccloLx_/x LicL-cty 


function of tiine following reactor shutdown. Figure 8 shows th< 
power to steady- state design power against time after reactor shutdown. The parameter 
in figure 8 is the reactor full-power operating time before shutdown. 



Figure 8. - Ratio of reactor decay power to steady-state poiver against time after shutdown with reactor 
fuii power operating time as a parameter. 


The lithium flow rates required to maintain a constant lithium temperature difference 
of 100° R (55.6 K) across the reactor core are shown against time after shutdown in fig- 
ure 9. These flow rates were calculated from equation (4), and are based on the reactor 
design operating power of 2. 17 megawatts-thermal. 

As explained in the section ANALYSIS, the upper limit of decay heat generation 
(which corresponds to 1 year of reactor full-power operation before shutdown) was used 
to determine the area requirements of the auxiliary heat exchanger and auxiliary radi- 
ator. From the procedure described in the section ANALYSIS, the following area and/or 
length requirements were determined: 

For the auxiliary heat exchanger. 

Length, Lj^^, ft (m) 1.81 (0.552) 

Heat-transfer surface area, A^ ft^ (m^) 0. 64 (0. 0594) 

(The cross-sectional configuration of the auxiliary heat exchanger is described in 
table II. ) 
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Figure 9. - Lithium flow rate against time after shutdown with reactor full power 
operating time as a parameter. Lithium temperature rise across core constant 
at 100° R (55. 5 K). 


For the auxiliary radiator, 

Planform length, ft (m) 6.0 (1.83) 

Effective radiation area (both sides), ft^ (m^) 37. 7 (3. 50) 

Planform width (tubes and fins), ft (m) 3. 14 (0. 958) 


(The cross-sectional configuration of the auxiliary radiator is described in table III. ) 
From the cross-sectional configuration and the value of the total mass of the 

auxiliary radiator (which includes tubes and fins only) was calculated to be approximately 
225 pounds (102 kg). 

The auxiliary radiator mass and the heat-transfer surface areas just listed were 
used to determine flow rates and temperatures in the auxiliary loop as a function of time 
after shutdown. Two sets of data were calculated. One set was based on a reactor full- 
power operating time of 1 year; the other set was for a full- power operating time of 
1 day. In the following paragraphs, we present data for each of these reactor operating 
times. 


Reactor Operating Time of 1 Year 

Calculated sodium flow rates and temperatures in the auxiliary loop are shown as a 
function of time after reactor shutdown in figures 10 to 14. These data (figs. 10 to 14) 
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Time after shutdown, sec 


Figure 10. - Total sodium flow rate in auxiliary loop and sodium flow rate in auxiliary 
heat exchanger against time after shutdown. Reactor operating time, lyear; reactor 
power, 2, 17 megawatts -thermal. 

are based on a reactor operating time of 1 year before shutdown. 

The total flow rate of sodium in the auxiliary loop and the flow rate which passes 
through the auxiliary heat exchanger are shown against time in figure 10. The difference 
between the two flow rates is equal to the auxiliary heat exchanger bypass flow rate; that 

^Na 6 ^Na 8 “ ^Na 7 ^ shutdown. 

At approximkely lOo’seconds after shutdown, the difference between the curves in 
figure 10 approaches zero. And the bypass flow rate at this time is essentially zero. 

The fact that the curves in figure 10 are tangent at some time after shutdown signifies 
that the area (or length) assigned to the auxiliary heat exchanger is indeed a minimum for 
the configuration assumed and described in table II. 
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Time after shutdown, sec 

Figure 11. ■- Totai sodium flow rate in auxiiiary ioop and sodium fiow rate in auxiliary 
radiator against time after shutdown. Reactor operating time, 1 year; reactor 
power, 2. 17 megawatts-thermal. 


The total sodium flow rate in the auxiliary loop and the flow rate in the auxiliary ra- 
diator are shown as a function of time after shutdown in figure 11. As indicated in fig- 
ure 11, for the time span from 0 to 166 seconds after shutdown, the total sodium flow 
passes through the auxiliary radiator . But for times greater than 166 seconds, a frac- 
tion of the total flow bypasses the auxiliary radiator . The bypass flow rate Wj^^ 2 is 
the difference between the total flow rate Wj^^^ g and the flow rate through the auxiliary 
radiator Wj^^ Figure 11 (b) shows that after 2 ^ hours (9000 sec), only about 20 per- 
cent of the total sodium flow goes through the auxiliary radiator . 

The mixed- mean temperatures of sodium at the auxiliary radiator inlet and auxiliary 
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Figure 12, - Mixed-mean temperature of sodium at auxiliary radiator inlet and auxiliary 
heat exchanger inlet against time after shutdown. Reactor operating time, 1 year; 
reactor power, 2. 17 megawatts -thermal. 


heat exchanger inlet are shown against time after shutdown in figure 12. The difference 
in temperature between these curves at any time t (i. e. , g - is constant and 

equal to 100° R (55. 6 K) throughout the shutdown period. The slopes of the curves in fig- 
ure 12 decrease with time and asymptotically approach a final value of zero. About 
166 seconds are required before the mixed-mean temperatures g and %a, 5^ 

reach their final values. For times greater than 166 seconds, the temperatures g 

and Tj^^ ^ are essentially constant and equal to 1760° R (978 K) and 1660° R (922 K), 
respectively. 

The shapes of the curves in figure 12 are influenced primarily by the effectiveness of 
the auxiliary radiator and the decay heat input to the auxiliary loop. The auxiliary radi- 
ator effectiveness increases with temperature and/or time after shutdown, whereas the 
decay heat generation rate (fig. 8) decreases with time. 

Figure 13 shows the sodium temperature at the auxiliary heat exchanger outlet 
against time after reactor shutdown. Shortly after a shutdown occurs, the sodium tem- 
perature at the auxiliary heat exchanger outlet is approximately 1300° R (722 K). It in- 
creases from this value at a significant rate during the first 400 seconds after shutdown; 
and thereafter it slowly approaches the maximum lithium temperature TLi,l in the pri- 
mary loop. 

During normal system operations (i. e. , prior to a shutdown), a small flow rate of 
sodium is required in the auxiliary heat exchanger in order to keep the sodium mixed- 
mean temperatures at the assumed standby level of 760° R (372 K). Even with the low 
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Time after shutdown, sec 

Figure 13. - Sodium temperature at auxiliary heat exchanger outlet against time after 
shutdown. Reactor operating time, 1 year; reactor power, 2, 17 megawatts-thermal. 

sodium flow rate which occurs during normal system operation, the liquid-sodium heat- 
transfer coefficient is relatively large. And as a result, the sodium temperature at the 
auxiliary heat exchanger outlet during normal system operation is near the maximum 
lithium temperature Tj^. j of 2200° R (1222 K). Hence, immediately after a shutdown, 
a sudden change must occur in the auxiliary heat exchanger outlet temperature. 

This momentary transient was not considered; and for this reason the temperature 
curve in figure 13 does not extend to the time zero point. 

Figure 14 shows the sodium temperature at the auxiliary radiator outlet as a func- 
tion of time after shutdown. The temperature Tj^^^ ^ reaches a maximum of 1660° R 
(922 K) at about 166 seconds after shutdown. (As shown in figure 11, this time also 
corresponds to the initiation of bypass flow around the auxiliary radiator. ) After about 
166 seconds, the temperature Tj^^^ ^ decreases with time. The reason for the gradual 
decrease in Tj,^^ ^ is that the sodium flow rate which goes through the auxiliary radi- 
ator gets smaller with time. As a result, the dwell time of the sodium in the auxiliary 
radiator gets larger with time after shutdown. 




Time after shutdown, sec 


Figure 14. - Sodium temperature at auxiliary radiator outlet against time after shutdown, 
Reactor operating time, 1 year; reactor power, 2. 17 megawatts-thermal. 


Reactor Operating Time of 1 Day 


The flow rates and temperatures of sodium in the auxiliary loop following a reactor 
full-power operating time of 1 day are shown in figures 15 to 19. Some of the data in 
these figures are similar to the data shown in figures 10 to 14 for a reactor operating 
time of 1 year. 

The total flow rate of sodium in the auxiliary loop and the flow rate through the 
auxiliary heat exchanger are shown against time in figure 15. A comparison of fig- 
ures 15 and 10 shows that the curves in figure 15 are not tangent at any point in time 
after shutdown; but in figure 10, the curves are tangent at approximately 106 seconds. 


The fact that the curves in figure 

2 

changer area (A. = 0. 64 ft or 


15 are not tangent signifies that the auxiliary hea 
2 

0. 0594 m ) is more than adequate to transfer the 
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Figure 15. - Total sodium flow rate in auxiliary loop and sodium flow rate in auxiliary 
heat exchanger against time after shutdown. Reactor operating time, 1 day; reactor 
power, 2. 17 megawatts-thermal. 


heat to the auxiliary loop for a reactor preshutdown operating time of 1 day. Conse- 
quently, a large fraction of the total sodium flow rate bypasses the auxiliary heat ex- 
changer throughout the shutdown. (The auxiliary heat exchanger bypass flow rate is the 
difference between the curves in fig. 15. ) 

Figure 16 shows the total sodium flow rate in the auxiliary loop and the sodium flow 
rate in the auxiliary radiator against time after shutdown. A comparison of figures 16 
and 11 shows that the bypass flow around the auxiliary radiator begins at approxi- 
mately the same time (about 166 sec after shutdown) for both reactor preshutdown oper- 
ating times. This indicates that the ratios of heat rejection rate to decay heat input rate 
for the two pre shutdown operating times are equal at about the same time after shutdown. 

Figure 17 shows the mixed-mean sodium temperatures at the auxiliary radiator 
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Time after shutdoivn, sec 

Figure 16, - Total sodium flow rate in auxiliary loop and sodium flow rate in auxiliary 
radiator against time after shutdown. fTeactor operating time, 1 day; reactor power, 
2. 17 megawatts -thermal. 



Time after shutdown, sec 

Figure 17. - Mixed-mean temperature of sodium at auxiliary radiator inlet and auxiliary 
heat exchanger inlet against time after shutdoivn. Reactor operating time, 1 day; 
reactor power, 2. 17 megavvatts-thermal. 


OK 



inlet and auxiliary heat exchanger inlet against time after shutdown. The temperatures 
in figure 17 increase with time and reach their final values at approximately 166 seconds 
after shutdown. After 166 seconds, the temperatures g and ^ are essen- 

tially constant and equal to 1710° and 1610° R (950 and 894 K), respectively. These final 
temperatures are approximately 50° R (27. 8 K) lower than those shown in figure 12. 

The sodium temperature at the auxiliary heat exchanger outlet is shown against time 
after shutdown in figure 18. This temperature curve is similar to that shown in figure 13. 

Figure 19 shows the sodium temperature at the auxiliary radiator outlet against time 
after shutdown. The temperature Tj^^ ^ reaches a maximum of 1610° R (894 K) at 
about 166 seconds after shutdown. (Fi^re 16 shows that this time also corresponds to 
the initiation of bypass flow around the axixiliary radiator. ) After 166 seconds, the tem- 
perature Tjg^ 2 decreases gradually with time. 

From a control standpoint, the range of sodium flow rates required during the shut- 
down period is important. The calculated data show that some of the sodium flow rates 
in the auxiliary loop change by two orders at magnitude (or more) during the first 9000 



Time after shutdown, sec 

Figure 18. - Sodium temperature at auxiliary heat exchanger outlet against time after 
shutdown. Reactor operating time, Iday; reactor power, 2. 17 megawatts-thermal. 
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Time after shutdown, sec 

Figure 19. - Sodium temperature at auxiliary radiator outlet against time after shutdown. 

Reactor operating time, 1 day; reactor power, 2. 17 megawatts-thermal. 

seconds after shutdown. To properly regulate these flow rates over such a wide range 
may require a special flow control valve arrangement. 

In the foregoing sections, we have shown how the flow rates and temperatures in the 
auxiliary loop vary with time after reactor shutdoAvn. These calculated data are based on 
a specified set of operating conditions along with selected cross-sectional configurations 
for the auxiliary heat exchanger and auxiliary radiator. The required areas of the auxil- 
iary heat exchanger and auxiliary radiator were determined to be relatively small. (The 

2 2 

auxiliary heat exchanger surface area was determined to be 0. 64 ft or 0. 0594 m ; 

2 2 

and the auxiliary radiator surface area was 37. 7 ft or 3. 50 m . ) 
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The reason that the auxiliary radiator is so small is due, in part, to the thermal 
capacitance of the auxiliary loop. Immediately following a shutdown, the decay heat gen- 
eration rates (fig. 8) are relatively large. And the thermal capacitance of the auxiliary 
loop is used during the early part of the shutdown period (i. e. , for about the first 60 sec 

or so) to absorb most of the decay heat. 


SUMMARY OF RESULTS 

An auxiliary radiator loop with sodium coolant was investigated as a possible means 
for removing decay heat from a proposed 2. 17-megawatt-thermal space power reactor. 
The' auxiliary loop considered for this study is made up of an auxiliary heat exchanger, an 
auxiliary radiator, a liquid-sodium pump, and the necessary flow control valves and pip- 
ing. The sodium-cooled auxiliary loop is connected to the reactor primary loop by the 
auxiliary heat exchanger. 

Based on a selected set of operating conditions and assumptions, minimum areas 
were determined for both the auxiliary heat exchanger and the auxiliary radiator. 

The minimum auxiliary heat exchanger surface area required to transfer the decay 
heat from the reactor primary loop to the axaxiliary loop was determined to be 0. 64 
square foot (0. 0594 m ). The auxiliary heat exchanger model used in this study is a 
simple counterflow type, consisting of two concentric cylindrical ducts. 

The minimum surface area required by the auxiliary radiator was determined to be 

2 

37. 7 square feet (3. 50 m ). Based on this surface area, the maximum radiator temper- 
ature was about 1760° R (978 K). The auxiliary radiator model consisted of 24 coplanar 
parallel tubes with central connecting fins of 100 percent effectiveness. The planform 
dimensions of the auxiliary radiator (assuming heat is rejected to space from both sides 
of the radiator) are: planform length, 6. 0 feet (1. 83 m); planform width, 3. 14 feet 
(0. 958 m). 

The configurations and minimum areas for the auxiliary heat exchanger and auxiliary 
radiator were used to determine sodium flow rates and temperatures throughout the 
auxiliary loop as a function of time after reactor shutdown. The sodium flow rates and 
temperatures were determined for preshutdown reactor operating times of 1 year and 
1 day. 

The calculated data show that some of the sodium flow rates in the auxiliary loop 
change by two orders of magnitude (or more) during the first 9000 seconds after shut- 
down. 
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A significant conclusion obtained from this preliminary study is that the size (and/ or 
mass) of the auxiliary loop is relatively small. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, May 26, 1970, 

120-27. 
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APPENDIX - SYMBOLS 


AA 

surface area of segment in ana- 
lytical model of auxiliary radi- 
ator, ft^; m^ 


thermal conductivity of auxiliary 
heat exchanger material (tanta- 
lum), Btu/(sec)(ft)(°R); 

\ 

surface area for heat transfer in 


kW/(m)(K) 


auxiliary heat exchanger based 
on inside surface of sodium an- 

Lrr 

Hx 

length of auxiliary heat exchanger, 
ft; m 


nular flow passage, ft^; 

^Rad 

length of auxiliary radiator, ft; m 


specific heat of radiator material 
(tantalum ) , Btu/ (Ibm ) (°R) ; 
kW-sec/(kg)(K) 

AM 

mass of radiator segment in ana- 
lytical model of auxiliary radi- 
ator, Ibm; kg 

S 

D„ 

specific heat, Btu/ (Ibm) (°R); 
kW-sec/(kg)(K) 

inside diameter of lithium flow 

Nu 

Pe 

Nusselt nxunber 
Peclet number 

a 

duct in auxiliary heat exchanger 

Pr 

Prandtl number 


(see fig. 4), ft; m 

Q 

reactor decay heat generation 


inside diameter of sodium annu- 


rate, Btu/sec; kW 

lar flow passage in auxiliary 

Re 

Reynolds number 


heat exchanger (see fig. 4), 
ft; m 

T 

T 

temperature, °R; K 
average coolant temperature in 

Dc 

outside diameter of sodium annu- 
lar flow passage in auxiliary 


radiator segment, °R; K 


heat exchanger (see fig. 4), 
ft; m 

Tr 

average material temperature of 
radiator segment, °R; K 

E 

effectiveness of counterflow auxil- 

t 

time after reactor shutdown, sec 


iary heat exchanger (eq. (8)) 

Ub 

overall heat-transfer coefficient 

h 

convective heat-transfer coeffi- 
cient in auxiliary heat ex- 
changer, Btu/(sec)(ft^)(°R); 

for auxiliary heat exchanger 
(eq. (3)), Btu/(sec)(ft2)(°R); 
kW/(m^)(K) 


kW/(m^)(K) 

w 

flowrate, Ibm/sec; kg/sec 


convective heat-transfer coeffi- 

e 

radiator surface emissivity 

cient in auxiliary radiator, 
Btu/(sec)(ft^)(°R); kW/(m^)(K) 

P 

density, Ibm/ft^; kg/m^ 
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a Stefan- Boltzmann constant, 

0. 476X10"^^ Btu/(sec)(ft^)(°R^); 
5.71x10'^^ kW/(m^)(K^) 

At auxiliary loop transport delay 

(eq. (15)), sec 

Subscripts: 

i inside surface of auxiliary radi- 

ator tubes (fig. 6) 

in inlet of segment in analytical 

model of auxiliary radiator 
(fig. 6) 

j segment number in analytical 

model of auxiliary radiator 

Li lithium 

Na sodium 

o outside surface of auxiliary radi- 

ator tubes (fig. 6) 


out outlet of segment in analytical 

model of auxiliary radiator 

(fig. 6 ) 

s sink 

1 reactor outlet (fig. 3) 

2 reactor inlet (fig. 3) 

3 auxiliary radiator bypass (fig. 3) 

4 auxiliary radiator outlet (fig. 3) 

5 auxiliary heat exchanger inlet 

(fig. 3) 

6 auxiliary heat exchanger bypass 

(fig. 3) 

7 auxiliary heat exchanger outlet 

(fig. 3) 

8 auxiliary radiator inlet (fig. 3) 
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